Kleinberg Navigation in Fractal Small World Networks 
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We study the Kleinberg problem of navigation in Small World networks when the underlying 
lattice is a fractal consisting of N 3> 1 nodes. Our extensive numerical simulations confirm the 
prediction that most efficient navigation is attained when the length r of long-range links is taken 
from the distribution P(r) ~ r~ a , where a = df, the fractal dimension of the underlying lattice. We 
find finite-size corrections to the exponent a, proportional to l/QnJV) 2 . 
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Recently Kleinberg has studied the problem of efficient 
navigation in Small World networks, based on local al- 
gorithms that rely on the underlying geography [l|, 
Consider for example a (i-dimensional hypercubic lattice, 
consisting of N nodes, where in addition to the links 
between nearest neighbors each node i is connected ran- 
domly to a node j with a probability proportional to r^ a , 
(here, and elsewhere, r,j = |r< — rj| denotes the Euclidean 
distance between nodes i and j). Suppose that a message 
is to be passed from a "source" node s to a "target" node 
t, along the links of the network, by a decentralized, or 
local algorithm (an algorithm that relics solely on data 
gathered from the immediate vicinity of the node that 
holds the message at each step), when the location of the 
target is publicly available. Kleinberg shows that when 
the exponent a = d an algorithm exists that requires less 
than (IniV) 2 steps to complete the task. If a ^ d, the re- 
quired number of steps grows as a power of N. Moreover, 
no local algorithm will do better, functionally, than the 
simple-minded greedy algorithm: pass the message for- 
ward to the neighbor node which is closest to the target 
(geographically) . 

Kleinberg observes Q that the above results generalize 
to "less structured graphs with analogous scaling prop- 
erties." Interest in such cases is practical, as the nodes 
of many real-life networks (routers of the Internet, pop- 
ulation in social nets, etc.) are not distributed regularly. 
Here we test this assertion for the case of fractal lattices, 
enhanced by the addition of random long-range links as 
in the original Kleinberg problem. We find that the re- 
sults indeed generalize to this case, and that most effi- 
cient navigation is achieved when the power exponent for 
the random connections is a = df, the fractal dimension 
of the underlying lattice. Our numerical analysis is sen- 
sitive enough to allow for a study of finite-size effects. 
For a lattice of N nodes optimal navigation is attained 
for an effective exponent a(N) that is smaller than the 
idealized limit of a = d{ (when N — > oo) by as much 
as l/(\nN) 2 . Thus, corrections are substantial even for 
very large lattices. 

Consider a fractal lattice, such as the Sierpinski car- 



pet [3j (Fig. [T]), where, in addition to the existing links, 
each node i is randomly connected to a single node j, 
selected from among all other nodes with probability 
Pij = r^ a j J2k r ik a ■ The sum in the denominator runs 
over all nodes k ^ i and is required for normalization. If 
the fractal is finite, consisting of N ^> 1 nodes, its linear 
size is L ~ N 1 ^ 1 . The normalization term then scales as 



k 



ik 



-a df 




The average distance between randomly chosen (s, im- 
pairs is ~ L, so in the absence of long contacts a message 
takes T ~ N 1 ^ 1 steps to be delivered Long range 
contacts reduce the 1/df exponent, but only when the 
exponent a — df does the expected delivery time scale 
slower than a power of N. The basic idea of Kleinberg's 
argument, applied to the case of fractals, is as follows [2]. 
For a = df, surround the target node t with m-shells of 



radii e r ' 



< r < e 



1,2,... Suppose that the 



message holder is currently in shell m, then the proba- 
bility that the node is connected by a long-range link to 
a node in shell m — 1, is, according to ([I]), 



1 



lnZ^=ffiX 



The probability to reach the next shell (m — 1) in more 
than x steps is then ip(x) = (1 — l/mL]^, and the average 
number of steps required to do so is 



(x) 



tp{x) dx ~ lnL 



'Electronic address: benavraham@clarkson.edu 



Since the largest shell is for e rn — L, the number of shells 
between the source and target is of the order of m — 
InL. Thus, the expected total number of steps required 
to reach the target is ~ (InL) 2 

For < a < df, surround the target node t by a 
ball of radius I — L 5 , (0 < <5 < 1). The probabil- 
ity that a randomly chosen node, i, has a long range 
contact to a site j within the ball, is, according to (JlJ, 
- r^ a jL dt - a < l/L dt '- a . Thus, the probability that 

i connects to any node in the ball does not exceed 
gd,f L dt- a = L Sdt-d {+a ^ Sincc s < 1; the source i ies al _ 

most surely outside of the ball (in the limit N — > oo). 
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FIG. 1: Small World net based on the Sierpinski carpet. 
Shown is a carpet of generation n = 3. The nodes (empty 
squares) are connected to their nearest neighbors (not shown) . 
In addition, each node i is connected to a random node j as 
described in the text. For the sake of clarity, only one such 
connection is shown as an example. 

Then, any £-step path between the source and target 
must contain at least one long range connection into a 
node inside the ball. But the probability that a node 
with such a connection is encountered within £ steps is 
smaller than I x 2/ df ~ df+Q . If this probability vanishes, 
as N — -> oo, the expected number of steps is bound to be 
at least I. This happens for 6 < (d{ — a) / (1 + d{) , so the 
expected number of steps exceeds L^ dt ~°^/^ 1+dl \ 

For a > G?f, the probability that a node has a long- 
range connection longer than L 1 , (0 < 7 < 1) scales, 
according to |T]), like 

/•oo -a 

J L -, a - d { 

Thus, the probability to jump a distance larger than L 1 
within L 13 steps, (0 < p < 1), is less than L/3 L7 (d f -a)^ 
If this probability vanishes, as A — > 00, then the total 
distance covered cannot exceed L^L 1 . Since the expected 
distance between source and target is of order L, we need 
j3 + 7 = 1. On the other hand, the probability for steps 
longer than L 1 vanishes when f3 + j(df — a) < 0. The two 
conditions yield (3 < (a — df)/(a — df+1), so the expected 
number of steps exceeds ~ L («-*)/(«-*+!). 

We have simulated navigation by the greedy algo- 
rithm in Small World nets based on the Sierpinski carpet 
(Fig. |lj) . A finite carpet, constructed recursively to gen- 
eration n, consists of A = 8™ nodes arranged within a 
square of side L — 3™. In each run random nodes are se- 
lected as the source (s) and target (t) and a path is sought 
between the two by the greedy algorithm. To minimize 
computer time and memory, each successive node and its 
random long-range connection is constructed only as it is 
reached by the message. In this way we were able to sim- 
ulate carpets of up to generation n = 19 {L w 1.2 x 10 9 , 



N w 1.4 x 10 17 ). For each value of the exponent a, the 
expected number of steps was obtained from averaging 
over 10,000 runs, apart from the largest lattice, of gen- 
eration n = 19, for which the number of runs per data 
point was reduced to 1,000. 

An example of the results obtained in this way is shown 
in Fig. [21 where we plot the logarithm of the "time" T 
- the average number of steps required by the greedy 
algorithm — as a function of a, for nets of generation 
n = 12, 15, and 18. Notice the parabolic shape of the 
curves, similar to what is observed for the regular square 
lattice [l[ , as well as the increase in T as the lattice grows 
larger. 
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FIG. 2: (Color online) Delivery time, T, as a function of the 
long-contact exponent, a. Plotted are results for the Sier- 
pinski carpet of generation n = 12, 15, and 18 (bottom to 
top). 

To extract further information, we have fitted second- 
order polynomials to the lnT(cv) curves of Fig. O The 
fits allow us to compute a(N)] the location of the mini- 
mum (or the optimal long-contact exponent) for a lattice 
of N nodes, as well as T(N); the minimal time required, 
on average, to deliver the message in a lattice of size 
N. In Fig. [3] we show the dependence of T(N) upon 
the lattice size. While it is quite hopeless to study the 
scaling T(N) ~ (mAf) 2 numerically (since it requires a 
double logarithmic function of N) , the slope of In T ver- 
sus In In N is seen to vary slowly and plausibly converging 
to 2, as predicted by Kleinberg. 

In Fig. [H we show the dependence of a(N) upon N. 
The optimal value of the exponent for long contacts varies 
very slowly with the size of the lattice. Inspired by the 
scaling T(N) ~ (In A^) 2 , we guess 

a(N) - a(oo) + A/ (hi A) 2 , (2) 

where A is a constant. This functional form, as well as 
the predicted limit of a(oo) = limAr^oo a(N) — df, for 
the case of an infinite lattice, is well supported by the 
data. The broadness of the lnT(a) curves (Fig. ^ makes 
it difficult to pinpoint their minima, a(A), leading to 
the large fluctuations evident in the plot. (On the other 
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FIG. 3: (Color online) Delivery time as a function of size. The 
logarithm of T is plotted against In In N. The data points are 
for carpets of generation n = 9, 10, . . . , 19 (from left to right). 
The straight line of slope 2, corresponding to the theoretical 
prediction T ~ (In TV) 2 , is shown for comparison. 



hand, T(N) can be determined quite accurately from the 
data, see Fig. [3]) 
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FIG. 4: (Color online) Optimal long-contact exponent. The 
exponent a(N), required for optimal navigation in a finite 
network of N nodes, is plotted against l/(lniV) 2 , for carpets 
of generation n = 9, 10, . . . , 19 (right to left). The arrow 
pointing at the expected limit of a(oo) = di = In 8/ In 3 is 
shown to guide the eye. 



In conclusion, we have studied navigation by the 
greedy algorithm on fractal Small World networks with 
random long-range connections taken from a power-law 
distribution. Our numerical results support the predic- 
tion of Kleinberg that optimal navigation occurs when 
the long-contact exponent is a = df, the fractal dimen- 
sion of the underlying lattice. 

An important effect are the corrections due to the fi- 
nite size of the networks involved, expressed in Eq. 
Even for our largest net, of N w 1.4 x 10 17 nodes, the 
correction to a is as big as 1.5%. More importantly, the 
correction is huge for sizes relevant to every-day life. For 
nets of 10 7 nodes (comparable to the size of the Inter- 
net, say), the correction is nearly 10% p}. It remains 
an open question whether corrections scale as l/(ln./V) 2 
generically, as observed in our case. 

We have also carried out simulations on a Small World 
network based on the Sierpinski gasket [3j. The re- 
sults are similar to the ones shown here for the car- 
pet, however, we observe a small but persistent dis- 
crepancy between the extrapolated a(co) ks 1.573 and 
d{ = In 3/ In 2 ks 1.585 of the Sierpinski gasket. For the 
ease of programming, we had embedded the gasket in a 
square lattice, thus distorting its original shape (of an 
equilateral triangle) to a right angled triangle. This in- 
troduces an anisotropy in the distribution of long-contact 
links: connections along the stretched direction are less 
favorable. We suspect that this anisotropy is the source 
of the discrepancy. Anisotropy effects will be considered 
in detail in future work. 
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